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. Abstract 
> , 

Q , A metliod for calculating exact propagators for those complex potentials with a real spec- 

' trum which are SUSY partners of real potentials is presented. It is illustrated by examples of 

. propagators for some complex SUSY partners of the harmonic oscillator and zero potentials. 



Recently a considerable attention has been paid to investigating different properties of non- 
J> I Hermitian Hamiltonians (see e.g. [1]). One of the reasons of that is an attempt to generalize the 



(N 



^ ' quantum mechanics by accepting non-Hermitian Hamiltonians with purely real spectrum (see e.g. 
fSJ ■ [1)2] and references therein). The propagator being the coordinate representation of the evolution 

^ . operator is one of the important objects in quantum mechanics since it permits to describe the 
\^ ! evolution of a quantum system for an arbitrary initial state. On the other hand the method of 
O I super symmetric quantum mechanics (SUSY QM) is one of the main methods for getting complex 
exactly solvable Hamiltonians [3, 4]. Moreover, because of its nice property to convert a non- 
diagonalizable Hamiltonian into diagonalizable forms and delete spectral singularities from the 



pH . continuous spectrum of a non-Hermitian Hamiltonian [5] it was presumed [6] that SUSY QM may 
^ ! become an essential ingredient of the infant complex quantum mechanics. The aim of this Letter is 
qh| to show that SUSY QM may be useful for finding propagators for those complex potentials which 
are SUSY partners of real potentials for which both the propagator and the Green's function are 
known thus giving an additional argument in favor of the above thesis. For simplicity we will 
consider only time-independent potentials although corresponding time-dependent technique is 
also available [7]. 

We consider the one-dimensional Schrodinger equation with a complex- valued potential Vc{x) 

[^^^-K]^x,t) = 0, (1) 

hc = -g^ + Vc{x), xeR. (2) 

Since the potential Vc{x) is time independent (stationary) solutions to equation (1) are expressed 
in terms of solutions of the stationary equation 

KipE{x) = EipE{x) (3) 

in the usual way, <l>(x, t) = e~'^^^(pE{x). But if initially the quantum system is prepared in the state 
(po{x) which is not a state with the definite value of the energy we need to know the propagator 



1 



for describing the evolution of a quantum system. We use the usual definition of the propagator 
Kc{x, y, t) as a solution to equation (1) with respect to variables x and t with the delta-like initial 
condition 

K,{x,y;0)^6{x-y). (4) 

If Kc{x, y, t) is given the function 

/oo 
K^{x,y;t)cpo{y)dy (5) 
-oo 

is a solution to equation (1) with the initial condition 0) = (fio{x). 

First of all we note that to be able to associate the function ^{x, t) with a state of a quantum 
system the integral in (5) should converge and both the function ipo{x) and ^{x,t) should belong 
to a certain class of functions. To present our method in its simplest form we will make several 
assumptions which shall simplify essentially our presentation keeping at the same time the essence 
of the method. Since in (5) the usual (Lebesgue) integration is involved it is natural to suppose 
that both (po{x) and t) are square integrable. This means that the Hilbert space, where the 
operator (non-Hermitian Hamiltonian) he associated with the differential expression —d^x/dx^ + 
Vc{x) 'lives', is the usual space L2(M) and the equation (3) creates an eigenvalue problem for 
he which is defined on a dense domain from L2(M). Similar eigenvalue problems were under an 
intensive study by mathematicians in the Soviet Union in the period between 50th and 70th of 
the previous century. Results of these investigations are mainly summarized in books [8, 9] to 
which we refer the interested reader where he can, in particular, find the strict definition of the 
spectrum, eigenfunctions, associated functions, domains of definition of operators created by non- 
Hermitian differential expressions and many other properties of differential equations and related 
non-selfadjoint operators. Here we would like to mention that the first essential result in this 
field was obtained by Keldysh [10] who proved the completeness of the set of eigenfunctions and 
associated functions for a non-self adjoint operator and results by Lidskiy [11] having the direct 
relation to the current paper. In particular, Lidskiy made a deep analysis of conditions on the 
potential Vc leading to an operator he which is uniquely defined by its closure and has a purely 
discrete spectrum with a complete set of eigenfunctions and associated functions. 

Our next essential assumptions is that he has a purely discrete spectrum, its set of associated 
functions is empty, it is diagonalizable, and its set of eigenfunctions n = 0, 1, . . . is complete 

in the space -L2(M). If h'^ is the adjoint differential expression it creates in the Hilbert space the 
adjoint operator with the eigenfunctions (()k{x) which also form a complete set in L2(M). Moreover, 
if En is an eigenvalue of he then E* (asterisk means the complex conjugation) is an eigenvalue of 
/i+ so that h'^(j)n = E*^(f)n- Note that neither {0^} nor {0„}, n = 0, 1, . . . form orthogonal systems 
but functions (pk are biorthogonal with 0„ and they can always be normalized such that (see e.g. 
[9]) 

/oo ^ 
(l)l{x)(l)n{x)dx = Snk ■ (6) 
-oo 

The completeness of the set of eigenfunctions of he means that any G L2(M) can be developed 
into the Fourier series over the set 0(a;) = X^^^o (^n<Pn{x)- Using the biorthonormality relation 

(6) we can find the coefficients in the usual way and put them back into the same relation thus 
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obtaining the symbolical form of the completeness condition of the set of eigenfunctions of he 



J2K{^)My) = s{x-y). (7) 



n=0 



Next we are assuming that the spectrum is real. Therefore the adjoint eigenvalue problem 
coincides with the complex conjugate form of equation (3) so that 4>n{x) — <Pn{^)- Under these 
assumptions equations (6) and (7) become (cf. with [12]) 

(l)n{x)(l)k{x)dx = 5nk : (8) 



= 5(x -y) . (9) 

n=0 

From here it follows the Fourier series expansion of the propagator in terms of the basis functions 

oo 

K,{x, y;t) = J2 M^)My)^~'''"' ■ (10) 

n=0 

Indeed, just like in the conventional Hermitian case the initial condition $(a;, 0) = (po{x) for 
function (5) with Kc of form (10) follows from (9) and the fact that ^{x,t) (5) satisfies equation 
(1) follows from the property of the functions 0„ to be eigenfunctions of he with the eigenvalues 

En- 

Especial role between all non-selfadjoint operators is played by pseudo-Hermitian operators first 
introduced by Dirac and Pauli and latter used by Lee and Wick [13] to overcome some difficulties 
related with using Hilbert spaces with an indefinite metric and their recent generalization (weak 
pseudo-Hermiticity) by Solombrino and Scolarici [14] since there are strict indications that these 
operators are the most appropriate candidates for replacing selfadjoint operators while generalizing 
the conventional quantum mechanics by accepting non- Hermitian operators [14, 15]. 

Another useful discussion is that the form (10) for the propagator may be interpreted as the 
coordinate representation of the abstract evolution operator.^ To show this we introduce ket- 
vectors (kets) \(pn) as eigenvectors of he and bra- vectors (bras) {4>n\ as functionals acting in the 
space of kets according to^ _ 

{(pn\(pk) ^ Snk ■ (11) 

Kets corresponding to the previous bras are just properly normalized eigenvectors of h'^ which 
is defined by the adjoint eigenvalue problem where Vc{x) is replaced by its complex conjugate 
V*{x) so that (11) is nothing but the same biorthogonality condition (8) written in the abstract 
representation. As usual the coordinate representation of the above abstract eigenvectors are 



"'^This property has been demonstrated in the report of the anonymous referee. 

^Without going into details we note that using the system one can construct the Hilbert space H so that 
the set of all finite linear combinations of is dense in H and formula (11) uniquely defines a functional in H, 
see e.g. [16]. 
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(/)n{x) = {x\(f)n) = {(f)n\x) where is an eigenvector of the coordinate operator. The completeness 
condition in the abstract form now reads 

oo 

J2\^n){M-^ (12) 
n=0 

and the formula 

oo 

K=J2\(l)n)En(^n\ (13) 

n=0 

presents the spectral decomposition of the Hamiltonian h^. Now in the known way [14, 15] one can 
introduce an automorphism rj to establish the property that he is (weakly) pseudo-Hermitian and 
construct a basis in which he takes a real form. We will not go into further details of the known 
properties of (weakly) pseudo-Hermitian operators since this is not the aim of this paper. The 
interested reader can consult papers [13, 14, 15] and the recent preprint [17] where biorthogonal 
systems are widely used in the study of different properties of non-Hermitian Hamiltonians. Our 
last comment here is that the abstract evolution operator given by its spectral decomposition 

oo 

m = |0n)e-^^"*(0n| (14) 

n=0 

written in the coordinate representation Kc{x,y,t) = {x\U(t)\y) is just the propagator (10). 

We would like to emphasis that conditions (8) and (9) have almost the usual form, only the 
complex conjugation is absent. Therefore they coincide with the corresponding equations for the 
Hermitian Hamiltonians in case when their eigenfunctions are real. 

The final assumption we make is that the Hamiltonian he is a SUSY partner of a Hermitian 
Hamiltonian ho with a purely discrete spectrum and a complete set of eigenfunctions ipn which 
always can be chosen real 

Q2 

ho = --g^ + Vo{x) , ho1pn{x) = En1pn{x) , ^/'n(x) = 7/'*(x) , 71 = 0,1,... 

SO that both the completeness and normalization conditions are given by equations (9) and (8) 
respectively with the replacement (pn ^ ipn- 

According to the general scheme of SUSY QM (sec e.g. [19]) operators h^ and he are (l-)SUSY 
partners if and only if there exists a first order differential operator L such that 

Lho^heL. (15) 

Operator L of the form 

L = -u'(x) /u{x) + d/d^ , (16) 

where the prime means the derivative with respect to x and the function u{x) is a solution to 
equation 

hou{x) = au{x) , (17) 
exists for any ho and a rather restricted set of operators he, the potentials of which are defined 

by 

Ve = Vo-2[\ogu{x)]". 
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The spectrum of /?.^. may either (i) coincide with the spectrum of or (11) may differ from it by 
one (real) level which is absent in the spectrum of ho. The case (i) may be realized only with a 
complex parameter a which is called the factorization constant. This statement is a consequence 
of a proposition proved in [4] . The case (ii) may be realized only for a real factorization constant 
since E = a is just the discrete level of he missing in the spectrum of ho and wc want that he has 
a real spectrum. Therefore in this case one has to choose u{x) as a linear combination of two real 
linearly independent solutions to equation (17). 

Together with operator L (16) we need also its 'transposed form' which we define as follows: 

L* = -u'(x) /u{x) - d/d^ . (18) 
Just like in the usual SUSY QM the following factorizations take place: 

L^L = ho-a, LL^ = he-a (19) 

which can easily be checked by the direct calculation. 

One of the main features of the method is that for the most physically interesting Hamiltonians 
ho operator (16) has the property L'^e(±cxo) = provided iPe{^oo). As a result the set of functions 

(t)n^NMn, n = 0,l,... (20) 

is complete in the space L^(R) in case (i). In case (ii) we have to add to this set the function 
= Na/u. The normalization coefficients Nn may be found by integration by parts in equation 
(8) and with the help of factorization property (19) which yields 

Nn^{En-a)-'/^. (21) 

The main result of the present Letter is given by the following 

Theorem 1. The propagator Kc{x,y;t) of the Schrodinger equation with the Hamiltonian he 

related with ho by a SUSY transformation is expressed in terms of the propagator Ko{x,y;t) of 

the same equation with the Hamiltonian ho and the Green's function Go{x, y; E) of the stationary 

equation with the same Hamiltonian as follows: 

in case (i) Kc{x,y,t) = KL{x,y,t) 

in case (ii) Kc{x,y,t) = KL{x,y,t) + 0a(ic)0„(y)e"*°* 

where KL{x,y,t) is the 'transformed' propagator 

/oo 
Ko{x,z,t)Go{z,y,a)dz . (22) 
-oo 

Here is defined by (16) and Ly is the same operator where x is replaced by y. 

Remark 1. We use the definition of the Green's function, which we denote Go{x,y, E) where E 
belongs to the resolvent set of ho, as a coordinate representation of the resolvent {ho — E)^^ (see 
e.g. [20]). The Green's function is a kernel of an integral operator acting in the space L2(M). 
There exists two equivalent representations for the Green's function: 
(a) In terms of the basis functions ipn{x) with eigenvalues E^, n = 0, 1, . . . 

G„(x,.,E) = f *fMf> (23) 

m=0 
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provided ipn{x) (ii"^ '"ea/ and 

(b) in terms of two special linearly independent solutions to the Schrddinger equation with the 
given value of E (see an example below). 

Proof. We prove the theorem only for case (ii) since case (i) does not have essentially new 
moments. 

Using equations (10), (20) and (21) we get 

K^{x, y, t) = L^Ly f2 "I^^^P^e-^^-' + 0.(x)0.(y)e-^"* . (24) 
Because of the equality 

/oo 
Ko{x,Z,t)iJrn{z)dz 
-oo 

equation (24) assumes the form 

K,(x, y, t) = L^Ly r Ko(x, z, t) M^)Myldz + 0a(x)0„(y)e— * . (25) 

m=0 ^ 

Here under the summation sign we recognize the Green's function (23) which ends the proof. □ 

We have to note that both in case (i) and in case (ii) a belongs to the resolvent set of Hq. Hence, 
the Green's function Go{x,y,a) entering in (22) and (25) is regular. Another useful comment is 
that almost the same proof is valid for the case when Hq (and hence he) has a continuous spectrum 
so that the theorem is valid for this case also provided he has no spectral singularities in the 
continuous part of the spectrum. 

To show how the theorem works we will consider a couple of typical examples. 

Consider first the harmonic oscillator Hamiltonian ho{x) = —d'^/dx^ + a;^/4 for which the 
propagator is well-known [21] 

1 i\(x'^+y'^) cos t — 2xy\ 

Ko{x,y,t) = e ^sint 

V 47r« sm t 

To apply our theorem wc need also the Green's function of the oscillator Hamiltonian at ii^ = a. 
To find it we use the definition of the Green's function in terms of two special solutions fi{x,E) 
and fr{x, E) of the equation 

(/lo - E)fi^r{.x, E) = Q, fi{-oo, E) = Q, /,(oo, ^) = 

which is 

G{x, y, E) = [fi{x, E)fr{y, E)Q{y - x) + fi{y, E)fr.{x, E)Q{x - y)]/W{fu fr) (26) 

where G is the Heaviside step function and W stands for the Wronskian. 
In the simplest case we can choose a — —1/2 and 

u{x) = e'^{C + eii{x/V2)), ImC^O 
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Functions // and fr from (26) at E — a — —1/2 read 

Mx, -1/2) = v^e^ (1 + erf {x/^/2)) , /.(x, -1/2) = v^e^ (1 - erf {x/V2)) . (27) 
The spectrum of the complex- valued transformed potential 

V,{x) = ^ - 1 + 2xQ^\x)e-'^ + 2gr'e""' , Qi{x) = ^[C + erf {x/^/2)] (28) 

consists of all oscillator energies £'„ = n+l/2,n = 0,l,... and one additional level -E_i/2 = —1/2 
with the eigenfunction 

0_i/2(x) = {2^:)-^l'^^C^u-\x) . (29) 
It is not difficult to check by the direct calculation that 

/oo 
(t>\j^{x)dx = l. (30) 
-oo 

Using Theorem 1 and equations (26) and (27) we obtain the propagator for the Hamiltonian 
with potential (28) 

K,{x, y, t) = -^^^^±ilL, j' Koix, z, ^e^ (1 + erf {z/V2))dz 

r Ko{x,z,t)e'^{l-eTf{z/V2))dz + cl>_y2{x)<f>-i/2{y)e''^^ . (31) 
V2u{y) Jy 

As the second example we consider a complex extension of the one-soliton potential 

Vcix) = . (32) 

cosh [ax + c) 

which is obtained from Vo{x) = with the choice 

u{x) — cosh(ax + c) , a — —a^ , Im c 7^ , Im a = 

as the transformation function. The Hamiltonian he with potential (32) has a single discrete level 
Eo — —a^ with the eigenfunction 



2 u{x) ■ 

Just like in the previous example this function is such that 



00 

2 



r_^2{x)dx^l. (33) 

To find the propagator Kc we are using the known propagator and Green's function for the free 
particle which are 



After some transformations similar to that described in details in [22] we get the propagator for 
the Hamiltonian he 



ae*" * 



K,{x,y;t) = ^==e « + [erf+ + erf_] (34) 



where the notation erf± = eri[a\/ii ± -^^{x — y)] is used and for c we take the value c = 

arctanh^^l^, lm(6) = 0. 

We have to note that formula (34) is in the prefect agreement with the result obtained by 
Jauslin [23] where the replacement t ^ it should be made since this author got the propagator for 
the heat equation with the one-soliton potential. Another comment we would like to make is that 
the Jauslin's method is based on an integral formula which relates solutions of two Schrdinger 
equations whose Hamiltonians are SUSY partners. Unfortunately, integrals which need to be 
calculated when the method is applied to the Schrodinger equation become divergent. There- 
fore the author found the propagator for the heat equation with the one-soliton potential. The 
Schrodinger equation may be considered as the heat equation with the imaginary time. In this 
respect the following question arises: whether the Jauslin's result after the replacement t ^ it 
gives the propagator for the Schrodinger equation with the one-soliton potential? We want to 
stress that the answer to this question is not trivial since such a replacement at the level of the 
Jauslin's integral transformation leads to divergent integrals and only the replacement in the final 
result gives the finite value for the propagator. Thus, our analysis shows that the answer to this 
question is positive. 

In conclusion we note that in this letter we presented a method for finding propagators for those 
non-Hermitian Hamiltonians with a purely real spectrum, SUSY partners of which have the known 
both the propagator and Green's function. The general theorem is illustrated by considering a 
complex anharmonic oscillator Hamiltonian and a complex extension of the one-soliton potential. 

The work is partially supported by the grants SS-5103.2006.2 and RFBR-06-02- 16719 and the 
Russian 'Dynasty' foundation (AMP). The authors would like to acknowledge that the form (14) 
of the evolution operator and its relation with the propagator was suggested by the anonymous 
referee. 
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